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Математическая модель прогнозирования влияния электромагнитного поля 
на устойчивость функционирования радиального подшипника, работающего 
на электропроводящей газовой смазке 


К. С. Ахвердиев, Е. О. Лагунова, К. А. Ванеев 
(Ростовский государственный университет путей сообщения) 


Представлен метод гидродинамического расчёта газодинамического радиального подшипника скольжения, 
работающем о на электропроводящей газовой смазке. Решается вопрос аналитического прогнозирования 
уст ОЙЧИВОСТИ работы подшипника скольжения. Дана оценка влияния числа Г. артмана, параметра сжимаемо- 
сти газа и параметра, характеризующего напряжённость электрического поля, на основные рабочие харак- 
теристики подшипника. 

Ключевые слова: газовая смазка, радиальный подшипник, автомодельное решение, сила трения, расход, 
Ггидродинамическое давление, несущая способность, Устойчивост! Ь. 


Введение. Как известно, подшипники с газовой смазкой в настоящее время широко используют- 
ся в машинах вращательного движения. Анализ существующих работ показывает, что в этой об- 
ласти есть ряд нерешённых проблем. Это, прежде всего, проблема, связанная с аналитическим 
прогнозированием устойчивости функционирования радиальных подшипников, работающих на 
электропроводящей газовой смазке при наличии магнитного поля. 

Цель настоящей работы — разработать аналитический метод расчёта радиального газо- 
динамического подшипника, работающего на электропроводящей газовой смазке. Дать оценку 
влияния электромагнитного поля на основные рабочие характеристики подшипника и на устойчи- 
вость его функционирования. 

Постановка задачи. Рассмотрим установившееся течение электропроводящей газовой смазки в 
зазоре радиального подшипника при наличии электромагнитного поля. Предполагаем, что под- 
шипник неподвижен, а шип вращается с постоянной угловой скоростью о. 

Будем исходить из уравнений «тонкого слоя» для вязкой несжимаемой жидкости при 
наличии электромагнитного поля. 

Основные уравнения и граничные условия. Запишем для нашей задачи уравнение движения 
газовой смазки, уравнение неразрывности и уравнение состояния [1]: 
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где и„,\ — компоненты вектора скорости; Р— гидродинамическое давление в смазочном слое; 
Е= [0,0,Е } — вектор напряжённости электрического поля; В = {В ‚0,0} — вектор магнитной 
индукции; о — электропроводимость газа; р’— плотность: А — удельная газовая постоянная; 


н — динамический коэффициент вязкости; /, 9 — полярные координаты. 
Функции Е (г',6) и В (г',6) считаем заданными, удовлетворяющими уравнениям Макс- 


велла [2, 3]: 
амВ = 0, гоёЁ = 0. (2) 
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Эти уравнения удовлетворяются при 


и г '(9 
Е = соп$Ё, В = и } (3) 





где Ф'(0)— заданная радиальная функция. 
В полярной системе координат с полюсом в центре шипа уравнения контура и подшипни- 
ка приближённо можно записать в виде 
=, Г'=Г+ес0$56, (4) 
где к — радиус шипа; и — радиус подшипника; е — эксцентриситет. 
Система уравнений (1) решается при следующих граничных условиях: 
и, =0, ц, =О при г' = г, +ес0$6; 
р (5) 
у, =0, и, = -ю^ приг’=/; Р(0)=Р(2п) = >". 
а 
Перейдём к безразмерным переменным по формулам: 


у. =юби, в =-@кц, г'= ПИ +бг, б=п-п,р(0)=^р(2п), 
Р=р.р, р'=р'р, р=р, Фф'=ф`ф. 


(6) 





Подставляя (6) в (1), с учётом (3) и (5) будем иметь с точностью до членов с[^) г с") : 


Го 





6? _ 1 ар д 
А в = 
ле №. ;(р)* ов (®')=0, б) 


и=0, ич=0 при 1+1с0$50; и =0, и=-р приг =0; =, 











2 52 
ОГ. № 
где Л = ы о — параметр сжимаемости газа; А = ор — параметр, характеризующий напря- 
р.б т 
ь 00?ф"? 
женность электрического поля, М = ая ЧИСЛО Гартмана. 
НА 


Точное автомодельное решение задачи. Точное автомодельное решение задачи (1)—(7) бу- 
дем искать в виде: 








оф оф Ф 
и=--—+И (г, 0), ру=—+И (г,0), Ф=——°—, 1(0)=1+псо0$6. 8 
й +0 р р т (6) ыы (8) 
Подставляя (8) в (7) получим следующие выражения: 
2 09, 2”, озу _ РР РА 
`гЗ г" 0? л @® 1(0) (9) 
0 ди 
= 
0г 00 
Для решения системы уравнений (9) воспользуемся р РИ условиями: 
И = (Е)пэтб; И =0(); Е= ; Ф=Ф(=), (10) 


56 й 
Подставляя (10) в (9), придём к следующей системе обыкновенных дифференциальных уравнений: 
рар _С, Ао, ай о 0, 

ло 1? № в Е С 
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азф ‚ ® = 

+ №ф. — =С., 
а РЕ (11) 
а?0 та 





Граничные условия для решения системы обыкновенных дифференциальных уравнений (11) бу- 
дут следующими: 


_Р (12) 


Решая линейные дифференциальные уравнения второго порядка с правой частью (11), с учётом 
граничных условий (12) найдём выражения для $" и 0: 


Ве яп (Ф5 УМЕ) [60$ (Ф.М) -1) Е 0$ (Ф УМЕ) С) 








2 2 $ 2 
Уп (ФУМ) № № № 
я п (ФУМ) [со (Фо УМ ) розм + 05 (ФУМ )С, -@) с05 (Ф МЕ) (№ +С) & 
0= и 
5 (Фо Ум ) № № № 
Из условия [06Е =0 найдём С,: 
ть А .р, 


№. [вов(о.ММ — 005 (в/м) +ап(.М) | 
в -2Мсоз(4,/М) +М№ + Моз(Ф, УМ т (Ф./М МФ, +9 (ФУМ) м 


Для нахождения безразмерного гидродинамического давления приходим к следующему 
уравнению: 








РР _ С.С, РА 











13 
ло Р'В ее 
Интегрируя уравнение (13) от 0 до 0, будем иметь 
я: 2^р, С С. ‚ АФ р 
= 1+ >< |@0. 14 
а Р. | г [|2 “з” Ш 


Решение интегрального уравнения (14) будем искать методом последовательных приближений, 
полагая в качестве первого приближения р» = 1. С точностью до С(Л?) для р; получим следующее 


выражение: 
и ее [4-8 С, а 


Из условия р(0) = р(2п) с точностью до О(п?) найдём С. : 





©. =-Аф, -А. 
Тогда выражение для гидродинамического давления р! с точностью до О(п?) примет вид: 
Ра 
р, = рр + р, ЛАпзт 0 + р,2ЛАф пт |. 
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Проекции главного вектора вещества смазок И/ = ню. + к, на шип определяется формулами: 
2п 2п 
К, =-, | Ряп@а6; К, = -^, | Рсо5@а@ 
0 0 


Заметим, что в силу периодичности давления Р и чётности с0$0 второй интеграл соответственно 
равен нулю. Первый интеграл проще вычислить по частям: 


т п ‚ГОР па 
К, =-1 [РА (-с056) = Рсоз6|" и, | —с05040=-и, | —<соз@ае. (15) 
0 0 [@5 0 [#5 


Подставляя выражение для гидродинамического давления в выражение (15), найдём выражение 
для несущей способности подшипника: 





р 
ШЕЕ, = = Г и (Ап+2Аф, + ЗАфоп). (16) 
а д 


Найдём силу трения по формуле: 











/ 12% в9=#® [оав- докл (А+ А) тм) 7 | (17) 
о ОГ бо б ФУМ [005 (Ф. УМ) +1) Р, 
Найдём расход смазки по формуле: 
О=р'® и] (=)4&= 
0 
©. 205 (ММ) М+М + оз( Ф.М) м &[№, чп (Ф.М) +9 (, УМ м) (18) 


= р’ гб 





р"озг,б 





№7 яп(Ф, /М < Ре яп(ф, /\ © 


Решение задачи об устойчивости. Решение этой задачи приводится на основе следующего 
безразмерного уравнения движения шипа: 


4? к [5\ аФ\ 
=- — с05Ф+& — |, 
ат? М \& ат 


4?Ф К, ‘(8 | 2% 92 
= этф-—| | — |. 
47? =5’Мг, Е\@ = ат дат 





(19) 








Здесь & = — относительный эксцентриситет; 7 = «6 Ё— время; и — угловая скорость вра- 
—Г 
1 0 
1/г 
щения шипа; е — эксцентриситет;  — радиус шипа; М — масса шипа; @& = |2 ; д — радиаль- 


ный зазор; 9 — ускорение силы тяжести; К, и А, — компоненты вектора поддерживающей силы; 
Ф — угол наложения. 

Явные выражения для А; и К„ приведённые выше, показывают, что А, и К, существенно 
зависят от следующих безразмерных параметров: параметра сжимаемости газа Л, параметра А, 
характеризующего напряжённость электрического поля, и параметра Л, характеризующего число 
Гартмана. Используя выражения для А, и А„ можно получить численное решение системы (19). 








4`=  а?’Ф ь 
Компоненты ускорения шипа ат? и т? представляют собой явные функции параметров 
Е @Ф ы 
= Ф, т!’ в, Л, А, М. Уравнения (19) записываются в стандартной форме первого порядка и 
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решаются с помощью многознакового метода, разработанного Гиром [4]. После получения реше- 
ь т Е @Ф е 

ния уравнений движения для заданных значений &, Ф, т’ в' Л, А, М устойчивость рас- 

сматриваемого движения определяется визуально по графику. Некоторые модели движения шипа 

приведены на рисунках 3 и 4. На рис. 3 представлено устойчивое движения шипа при следующих 


условиях п = 0,01; =, =0,5; 5= 2,4; Ф=0,001; Л=0,25; А=0,6; М =0,4. На рис. 4 представле- 





ь о 
но неустойчивое движение шипа (=, =0,8; —=2,5; Ф=0,001; л=0,25; А=0,25; А=0,1; 
[п 
№ =0,2). 
0,4 0,5 05 04 
0,6 0,3 09 10 09 0,8 я 0,6 Е 0,6 о а 


А 


Рис. 1. Зависимость безразмерной несущей способности . 
Рис. 2. Зависимость безразмерной силы трения /[ от числа 
от числа Гартмана, параметра сжимаемости газа и Гартмана и параметра А, характеризующего напряжённость 


электрического поля 





К а 
параметра А, характеризующего напряжённость электриче- 
ского поля. 1 — Л = 0,25; 2 —Л = 0,5 
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Относительное смещение из 
положения равновесия 





Относительное смещение из области равновесия 
> 


1 2 3 4 5 6 7 8 3 
Безразмерное время, Т 


123465 67 8 9 
Безразмерное время, Т 


Рис. 3. Движение центра шипа при моделировании устой- Рис. 4. Движение центра шипа при моделировании неустой- 
чивого СОСТОЯНИЯ чивого СОСТОЯНИЯ 
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1 234065 6 7 
1—Л= 0,0001, А=0, М= 0; 2 —л= 0,0007, А=0, М=0; 3 —1= 0,01, А=0, М= 0; 4 —п= 0,0001, А = 0,5, №М= 0,4 
Рис. 5. Область устойчивости движения шипа 
Результаты для некоторых значений Ф, Е, Л, Аи № приведены на рис. 5. Все точки, ко- 
торые лежат ниже кривых, соответствуют устойчивому движению, а все точки, которые лежат 
выше этих кривых, соответственно — неустойчивому движению шипа. Графики устойчивости 
определяются такими точками ( в)/6 ), которые являются неграничными, т. е. соответствуют орби- 


те шипа. 
Выводы. Результаты численного анализа аналитических выражений (16)—(18) для основных ра- 
бочих характеристик радиального подшипника показывают: 

— Число Гартмана не оказывает существенного влияния на несущую способность подшип- 
ника (рис. 1). С увеличением параметра А, характеризующего напряжённость электрического по- 
ля, и параметра сжимаемости газа имеет место значительное увеличение несущей способности. 

— С увеличением параметра А, характеризующего напряжённость электрического поля, 
безразмерная сила трения возрастает (рис. 2), а с увеличением числа Гартмана уменьшается. 
Принимая во внимание решение для малого числа приближений, параметр сжимаемости газа не 
оказывает существенного влияния на безразмерную силу трения. 

— Сувеличением числа Гартмана расход смазки резко уменьшается. При малых значениях 
с увеличением параметра А, характеризующего напряжённость электрического поля, расход смаз- 
ки незначительно возрастает. Принимая во внимание решение для малого числа приближений, па- 
раметр сжимаемости газа не оказывает существенного влияния на расход смазочного вещества. 

— Из приведённых на рис. 5 зависимостей следует, что в случае электропроводящей газо- 
вой смазки область устойчивости движения шипа расширяется. 
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